Defining window-boundaries for genomic analyses using smoothing spline techniques

Background

High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome.

Results

Simulations applying this method were performed to identify selection signatures from pooled sequencing F_ST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing F_ST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach.

Conclusions

We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on F_ST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.     

A GWAS assessment of the contribution of genomic imprinting to the variation of body mass index in mice

Background

Genomic imprinting is an epigenetic mechanism that can lead to differential gene expression depending on the parent-of-origin of a received allele. While most studies on imprinting address its underlying molecular mechanisms or attempt at discovering genomic regions that might be subject to imprinting, few have focused on the amount of phenotypic variation contributed by such epigenetic process. In this report, we give a brief review of a one-locus imprinting model in a quantitative genetics framework, and provide a decomposition of the genetic variance according to this model. Analytical deductions from the proposed imprinting model indicated a non-negligible contribution of imprinting to genetic variation of complex traits. Also, we performed a whole-genome scan analysis on mouse body mass index (BMI) aiming at revealing potential consequences when existing imprinting effects are ignored in genetic analysis.

Results

10,021 SNP markers were used to perform a whole-genome single marker regression on mouse BMI using an additive and an imprinting model. Markers significant for imprinting indicated that BMI is subject to imprinting. Marked variance changed from 1.218 ×10⁻⁴ to 1.842 ×10⁻⁴ when imprinting was considered in the analysis, implying that one third of marked variance would be lost if existing imprinting effects were not accounted for. When both marker and pedigree information were used, estimated heritability increased from 0.176 to 0.195 when imprinting was considered.

Conclusions

When a complex trait is subject to imprinting, using an additive model that ignores this phenomenon may result in an underestimate of additive variability, potentially leading to wrong inferences about the underlying genetic architecture of that trait. This could be a possible factor explaining part of the missing heritability commonly observed in genome-wide association studies (GWAS).     

Genome-wide differentially methylated genes in prostate cancer tissues from African-American and Caucasian men

Increasing evidence suggests that aberrant DNA methylation changes may contribute to prostate cancer (PCa) ethnic disparity. To comprehensively identify DNA methylation alterations in PCa disparity, we used the Illumina 450K methylation platform to interrogate the methylation status of 485,577 CpG sites focusing on gene-associated regions of the human genome. Genomic DNA from African-American (AA; 7 normal and 3 cancers) and Caucasian (Cau; 8 normal and 3 cancers) was used in the analysis. Hierarchical clustering analysis identified probe-sets unique to AA and Cau samples, as well as common to both. We selected 25 promoter-associated novel CpG sites most differentially methylated by race (fold change > 1.5-fold; adjusted P < 0.05) and compared the β-value of these sites provided by the Illumina, Inc. array with quantitative methylation obtained by pyrosequencing in 7 prostate cell lines. We found very good concordance of the methylation levels between β-value and pyrosequencing. Gene expression analysis using qRT-PCR in a subset of 8 genes after treatment with 5-aza-2′-deoxycytidine and/or trichostatin showed up-regulation of gene expression in PCa cells. Quantitative analysis of 4 genes, SNRPN, SHANK2, MST1R, and ABCG5, in matched normal and PCa tissues derived from AA and Cau PCa patients demonstrated differential promoter methylation and concomitant differences in mRNA expression in prostate tissues from AA vs. Cau. Regression analysis in normal and PCa tissues as a function of race showed significantly higher methylation prevalence for SNRPN (P = 0.012), MST1R (P = 0.038), and ABCG5 (P < 0.0002) for AA vs. Cau samples. We selected the ABCG5 and SNRPN genes and verified their biological functions by Western blot analysis and siRNA gene knockout effects on cell proliferation and invasion in 4 PCa cell lines (2 AA and 2 Cau patients-derived lines). Knockdown of either ABCG5 or SNRPN resulted in a significant decrease in both invasion and proliferation in Cau PCa cell lines but we did not observe these remarkable loss-of-function effects in AA PCa cell lines. Our study demonstrates how differential genome-wide DNA methylation levels influence gene expression and biological functions in AA and Cau PCa.     

Stimulation of allergen-loaded macrophages by TLR9-ligand potentiates IL-10-mediated suppression of allergic airway inflammation in mice

BackgroundPreviously, we demonstrated that OVA-loaded macrophages (OVA-Mφ) partially suppress OVA-induced airway manifestations of asthma in BALB/c mice. In vitro studies showed that OVA-Mφ start to produce IL-10 upon interaction with allergen-specific T cells, which might mediate their immunosuppressive effects. Herein, we examined whether IL-10 is essential for the immunosuppressive effects of OVA-Mφ in vivo, and whether ex vivo stimulation of the IL-10 production by OVA-Mφ could enhance these effects.MethodsPeritoneal Mφ were loaded with OVA and stimulated with LPS or immunostimulatory sequence oligodeoxynucleotide (ISS-ODN) in vitro. The increase of IL-10 production was examined and, subsequently, ex vivo stimulated OVA-Mφ were used to treat (i.v.) OVA-sensitized mice. To further explore whether Mφ-derived IL-10 mediates the immunosuppressive effects, Mφ isolated from IL-10^-/- mice were used for treatment.ResultsWe found that stimulation with LPS or ISS-ODN highly increased the IL-10 production by OVA-Mφ (2.5-fold and 4.5-fold increase, respectively). ISS-ODN stimulation of OVA-Mφ significantly potentiated the suppressive effects on allergic airway inflammation. Compared to sham-treatment, ISS-ODN-stimulated OVA-Mφ suppressed the airway eosinophilia by 85% (vs. 30% by unstimulated OVA-Mφ), IL-5 levels in bronchoalveolar lavage fluid by 80% (vs. 50%) and serum OVA-specific IgE levels by 60% (vs. 30%). Importantly, IL-10^-/-Mφ that were loaded with OVA and stimulated with ISS-ODN ex vivo, failed to suppress OVA-induced airway inflammation.ConclusionsThese results demonstrate that Mφ-derived IL-10 mediates anti-inflammatory responses in a mouse model of allergic asthma, which both can be potentiated by stimulation with ISS-ODN.     

Comparison of classification methods for detecting associations between SNPs and chick mortality

Multi-category classification methods were used to detect SNP-mortality associations in broilers. The objective was to select a subset of whole genome SNPs associated with chick mortality. This was done by categorizing mortality rates and using a filter-wrapper feature selection procedure in each of the classification methods evaluated. Different numbers of categories (2, 3, 4, 5 and 10) and three classification algorithms (naïve Bayes classifiers, Bayesian networks and neural networks) were compared, using early and late chick mortality rates in low and high hygiene environments. Evaluation of SNPs selected by each classification method was done by predicted residual sum of squares and a significance test-related metric. A naïve Bayes classifier, coupled with discretization into two or three categories generated the SNP subset with greatest predictive ability. Further, an alternative categorization scheme, which used only two extreme portions of the empirical distribution of mortality rates, was considered. This scheme selected SNPs with greater predictive ability than those chosen by the methods described previously. Use of extreme samples seems to enhance the ability of feature selection procedures to select influential SNPs in genetic association studies.     

Enterococcal surface protein Esp is not essential for cell adhesion and intestinal colonization of Enterococcus faecium in mice

Background  

Enterococcus faecium has globally emerged as a cause of hospital-acquired infections with high colonization rates in hospitalized patients. The enterococcal surface protein Esp, identified as a potential virulence factor, is specifically linked to nosocomial clonal lineages that are genetically distinct from indigenous E. faecium strains. To investigate whether Esp facilitates bacterial adherence and intestinal colonization of E. faecium, we used human colorectal adenocarcinoma cells (Caco-2 cells) and an experimental colonization model in mice.     

A novel liver specific isoform of the rat LAR transcript is expressed as a truncated isoform encoded from a 5'UTR located within intron 11

Background  

The leukocyte common antigen related receptor (LAR) protein has been shown to modulate the signal transduction of a number of different growth factors, including insulin and insulin-like growth factor 1. Splice variants exhibit differing roles and are expressed according to tissue type and developmental stage.     

Removing celiac disease-related gluten proteins from bread wheat while retaining technological properties: a study with Chinese Spring deletion lines

Background  

Gluten proteins can induce celiac disease (CD) in genetically susceptible individuals. In CD patients gluten-derived peptides are presented to the immune system, which leads to a CD4⁺ T-cell mediated immune response and inflammation of the small intestine. However, not all gluten proteins contain T-cell stimulatory epitopes. Gluten proteins are encoded by multigene loci present on chromosomes 1 and 6 of the three different genomes of hexaploid bread wheat (Triticum aestivum) (AABBDD).